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ABSTRACT 

This work combines new observations of NIR starlight hnear polarimetry with previously simulated 
observations in order to constrain dynamo models of the Galactic magnetic field. Polarimetric obser- 
vations were obtained with the Mimir instrument on the Perkins Telescope in Flagstaff, AZ, along a 
line of constant Galactic longitude {i = 150°) with 17 pointings of the 10' x 10' field of view between 
—75° < 6 < 10°, with more frequent pointings towards the Galactic midplane. A total of 10,962 
stars were photometrically measured and 1,116 had usable polarizations. The observed distribution 
of polarization position angles with Galactic latitude and the cumulative distribution function of the 
measured polarizations are compared to predicted values. While the predictions lack the effects of tur- 
bulence and are therefore idealized, this comparison allows significant rejection of AO-type magnetic 
field models. SO and disk-even halo-odd magnetic field geometries are also rejected by the observa- 
tions, but at lower significance. New predictions of spiral-type, axisymmetric magnetic fields, when 
combined with these new NIR observations, constrain the Galactic magnetic field spiral pitch angle 
to -6° ±2°. 

Subject headings: ISM: dust - ISM: magnetic fields - Polarization - Radiative transfer 



1. INTRODUCTION 

The Galactic magnetic field is an important aspect of 
the interstellar medium, but even its large-scale structure 
is still uncertain. What symmetry properties does the 
field exhibit? Is the magnetic field axisymmetric, or are 
there field direction reversals at some radii? How is the 
magnetic field sustained? As a first step to addressing 
these questions, this work compares new near-infrared 
(NIR) starlight polarimetry to predictions of the polar- 
ization properties for different magnetic field geometries 
(|Pavelll2011l ). Through strong observational constraints, 
the mechanism that sustains the Galactic magnetic field 
can also be constrained. In addition to the large-scale 
symmetry of the Galactic magnetic field, these observa- 
tions are used to constrain the spiral-type magnetic pitch 
angle of the Galaxy. 

Magnetic fields remain one of the most difficult aspects 
of the nearby universe to probe. Every method for mea- 
suring the interstellar magnetic field requires some spe- 
cial circumstance (e.g., a background polarized source 
for Faraday rotation, high column density for polarized 
thermal emission or Zeeman splitting, a background star 
for interstellar polarization of background starlight). To 
date, the study of the Galactic magnetic field has been 
dominated by radio wavelength Faraday rotation stud- 
ies. However, by combining many different techniques, 
observational constraints can be placed on the structure 
of the Galactic magnetic field which will drive, for exam- 
ple, dynamo theories forward. 

Over the last few decades, progress has been 
made in understanding the Galactic magnetic field, 
throu gh, especially, improvements in numerical s i mula- 
tions (Elstner et al. 1992VBrandenburg et al."1992','1993; 
Ferriere & Schmitt 2000 ; Kleeorin g t al.. 2 002, 200i 
Moss fc Sokolofa 120081: iHanasz'et al.ri2009l: iMoss et alj 
20101 ). However, while our understanding of fundamental 
dynamo physics grew, there were few observational pre- 

|pavelmi@bu.edu; clemens@bu.edu; apinnick@bu.edu| 



dictions. Observational constraints of large-scale mag- 
netic field models have almost exclusively relied on Fara- 
day rotation of emission from Galactic pul s ars and po- 
larized extragalactic sources (iSmithl 119681: iManchesteil 
[1971 IRand fc Kulka^ 119891: IMen et all I2008D which 
only probe the line-of-sight component of the magnetic 
field. .Heile s (1996) used optical polarimetry to con- 
strain the Galactic magnetic pitch angle (defined as 
p = tan~^[Br I B^) for a spiral-like magnetic field pat- 
tern), but, beyond this, starlight polarization has not 
been quantitatively employed as a tool for constraining 
the large-scale Galactic magnetic field structure. 

Polarization of background starlight has bee n an as- 
trophysical tool s i nce th e pioneering works of iHiltneii 
(l94cj) and iHalll (|19490 . and early theoretical work 
linked this optical polarization with interstellar mag- 
netic fields (Davis & Grccnstcin 1951). As currently un- 
derstood (Lazarian 2007), interstellar polarization arises 
from dichroic extinction of unpolarized starlight passing 
through regions of elongated dust grains which are direc- 
tionally aligned with their long axes preferentially per- 
pendicular to the local magnetic field. Assuming back- 
ground stars emit all polarizations equally (i.e., unpolar- 
ized), photons polarized parallel to a grain's long axis see 
the largest grain cross section and are preferentially ex- 
tincted. The result is linear polarization whose stronger 
electric field direction is parallel to the direction of the 
magnetic field, as projected on the sky. 

Much of the previous starlight polarization work has 
been done at optical wavelengths w here t he polarization 
signal is strongest (Serkowsk i et al.lll975t ). However, op- 
tical wavelengths suffer appreciable dust extinction, es- 
pecially through the Galactic disk, and typically can- 
not p robe stars farther than 2-3 kpc (Fosalba 2002| iHanl 
120081 ). At longer wavelengths, starlight is less extincted 
by dust and polarizations can be measured to stars far 
beyond optical limits. 

In addition to extinction, nearby magnetic structures 
can dominate the observed optically-traced morphology. 
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Loop I is a supernovae remnant of 116° diameter on 
the sky, approximately 130 pc from the S un, centered 
at ^ = 329°, b = 17.5° (jBerkhuiisenlflOTl . and promi- 
nently seen in the optical polarization compilation of 
iHeilesI ()200Cll ). This nearby object obscures the signature 
from the large-scale magnetic field structure. Faraday 
rotation studies of the Galactic magnetic field may also 
be affecte d by similar foreground objects (e.g., supernova 
remnants: IWolleben et"alT[20Tot) . By observing in a di- 
rection without any obvious disturbances, the quiescent 
Galactic magnetic field may be reliably probed by NIR 
starlight polarimetry to distances of several kiloparsecs 



Table 1 

Observations and Starlight Polarization Detections 



(IClemens et al.l l2012aD. 
IPavell ^ 20111) used existing dynamo-driven Galactic 



magnetic field models and empirical dust distributions 
to predict the NIR observational signatures of different 
magnetic field models: SO (even), AO (odd), disk-even 
halo-odd (DEHO), and simple analytic axisymmetric. 
From these predictions, the shape of the curve of po- 
larization Galactic position angle (hereafter GPAQ) with 
Galactic latitude (6) and cumulative distribution func- 
tions (CDFs) of the normalized degree of polarization (P) 
were proposed as tools for testing predictions of the large- 
scale structure of the Galactic magnetic field against ob- 
servations. ^^^^ 

Here, the predictions of lPavell (|2011l ) are tested against 
observations of NIR starlight polarimetry with the goal 
of constraining possible large-scale magnetic field geome- 
tries. In §2, the observations and data reduction are de- 
scribed. A summary of the simulated polarization mea- 
surements from existing Galactic dynamo and dust mod- 
els is presented in §3, along with a comparison between 
the predicted and observed polarimetric properties. The 
results of this analysis are discussed in §4, and conclu- 
sions are presented in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations were made with the Mimir instrument 
(jClemens et al.|[2007l) in H-band (1.6^m) linear imaging 
polarimetry mode on the 1.8m Perkins Telescope outside 
Flagstaff, AZ, on several nights from 2007 through 2009. 
Mimir uses a cold, stepping half-wave plate (HWP) and 
a cold, fixed wire grid to analyze starlight polarization 
across a 10' x 10' field of view. The detector is a 1024 x 
1024 InSb Aladdin III array. 

The observations were taken along a line of constant 
Galactic longitude, i = 150°, for -75° < b < 10° in 
steps of Ab = 5° to enable comparison with the IPavell 
(|201lD predictions. Observations taken at & = —5° were 
not used because the bright star HD 23049 contaminated 
the images, instead the neighboring b = —6° field was 
substituted. Galactic latitudes 6 = 0°, —35°, —50°, and 
—55° were not observed because of poor weather and the 
Perkins telescope is not able to observe North of 6 = 10° 
or South of = —75° at this Galactic longitude due 
to telescope mount limits. Supplemental measurements 
were obtained near the Galactic midplane at & = ±1.25° 
and ±2.5° to provide higher latitude resolution there. 

For each pointing, the HWP is rotated to 16 differ- 
ent position angles, equivalent to five measurements at 
an instrument position angle (IPA) of 0° and four IPA 
measurements each at 45°, 90°, and 135°. This is done 

^ Throughout this work, polarization position angles will be mea- 
sured in the Galactic coordinate system, measured East {+£) from 
Galactic North (+6). 
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Significant detections defined as P/a-p > 3 
^ Marginal detections defined as 1.25 < P/up < 3. 
° Upper limits defined as P/crp < 1.25 and and ap < 1.0%. 

toward six different sky dither positions, following a ro- 
tated hexagon pattern on the sky, for a total of 102 im- 
ages per pointing. Each target field consisted of four sets 
of observations with 10s integrations per HWP position, 
for a total of 68 minutes of integration time per field. In- 
tegration times were set to ensure that the polarimetric 
uncertainty at < 14 was less than 1%. The dates of 
each Galactic latitude observation are shown in column 
(2) of Table [TJ Most of the stars foun d have comple- 
ment ary NIR photometry from 2MASS (jSkrutskie et al.l 
[2006h . 

Polarization standards were also observed for calibra- 
tion. For the 2007 observations, the unpolarized star HD 
42807 (Gchrcls 1974) and polarized standard star HD 
30675 ( Whittct , 199211 were observed for calibration. For 
all other observations, the Galactic globular cluster NGC 
5466 was observed as a collection of unpolarized stars and 
the polarization st andard stars Elias 14, 22, and 25 in 
the Scorpius field dWhittetl [I992I: IClemen s et al. 2012b[) 
were observed. Detector linearity was characterized by 
observing the constant surface brightness of an in-dome 
screen for a range of exposure times, and each of the 16 
HWP position angles had a unique dome fla t measured, 
described in detail in 'Cleme ns et al.l ()2012al ) . 

Data reduction was performed with two custom IDL 
software packages (available through the Mimir web- 
sitj3), whose details can be found in IClemens et al.l 
(|2012aD . The first software package performs quality 
testing and applies linearity, dark, and polarimetric flat 
field corrections. The second package calculates astro- 
metric solutions for each image, coadds images, extracts 
photometry, calculates H-band starlight polarizations, 
and applies instrumental polarization corrections. Since 
each of the fields was observed with four pointings at 
six dither positions per pointing in sixteen unique HWP 
positions per dither, a total of 24 ten-second exposures 
were taken through each HWP position of each field (the 
HWP IPA= 0° position had 48-ten second exposures for 

^ http:/ /people. bu.edu/clemens/mimir/software. html 
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each field). For each field, these 24 exposures were coad- 
ded to form sixteen master HWP im ages per field. P SF- 
assisted aperture photometry (.Clemens et al.l[2012al ) of 
stars in these master HWP images were used to calculate 
the polarization properties of each star. The final output 
is a catalog of stellar polarizations (polarization percent- 
age, equatorial position angle. Stokes U, Stokes Q, and 
associated uncertainties) with Mimir H-band photome- 
try and 2MASS JHKg photometry, where available. The 
sixteen master HWP images were also coadded to create 
deep photometric images for each field. 
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Figure 1. Mimir coadded H-band intensity image of the 10' X 10' 
field toward b = —6° with overlaid significant polarization vectors. 
The length of each vector represents the degree of polarization 
relative to the 2% scale in the lower left corner, and the vector 
orientation shows the equatorial polarization position angle. The 
diagonal light gray line shows a line of constant Galactic latitude 
at 6 = -6° . 

Figure [T] shows the coadded H-band intensity image 
obtained toward b = —6°, with measured polarizations 
shown as vectors. For detected polarizations, the length 
of the vector represents the degree of starlight polariza- 
tion and its orientation represents the equatorial position 
angle of its electric field vector. A line of constant Galac- 
tic latitude is shown for reference. 

The £ = 150° region of sky was chosen for two rea- 
sons. First, the outer Galaxy provides a relatively qui- 
escent region (compared to the inner Galaxy) where the 
large-scale magnetic field can be observed. The inner 
Galaxy hosts much of the Galaxy's star formation and ex- 
hibits extensive overlapping supernova remnants, which 
can distort the magnetic field on parsec to kiloparsec 
scales. Second, the i = 150° field s overlap completed 
SDSS SEGUE (lYannv et al.l 120091) imaging and spec- 
troscopy fields. Subsequent analysis of the photometry 
and spectroscopy obtained by SEGUE may be used to 
inform distance estimates for the stars observed here to 
enable decomposing the magnetic field projections with 
distance, though this is beyond the scope of this current 
work. 

A total of 10,962 stars were detected in the com- 
bined photometric images, of which 332 had significant 
{P/crp > 3) polarization detections, 784 had "marginal" 
(1.25 < P/(Jp < 3) detections, and 194 had signifi- 
cant upper limits {P/up < 1.25 and ap < 1%). All of 
these 1,310 polarization 'targets' are listed in Table [5] by 
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Figure 2. Distribution of polarization uncertainty as a function of 
H-band magnitude. The solid line traces the median polarization 
uncertainty with apparent magnitude. The median polarimetric 
uncertainty at = 14 is shown by the dashed lines. 

Galactic coordinates, degrees of polarization, GPAs, and 
2MASS H- and K-band magnitudes. The lowest, median, 
and highest significant detected polarizations are 0.19%, 
1.51%, and 12.13%, respectively. All polarization target 
coordinates were passed through SIMBAD to search for 
previous spectral classifications, but only eight out of the 
1,310 polarization targets had entries. Of these matches, 
three stars had spectral classifications (A2V, A3V, and 
M3V) and one was ide ntified as the X -ray source IRXS 
J042738.8-1-555831 (Vo ges et al.llT999; ). while the other 
four were only identified as stars. The SIMBAD search 
result suggests that this is a quiescent direction in the 
Galaxy suitable for probing the large-scale Galactic mag- 
netic field. 

To aid in determining the completeness of these obser- 
vations, the distribution of H-band magnitudes for all de- 
tected stars was examined. From these distributions, the 
photometric observations appear complete to H ^ 16.5 
mag and the polarimetric observations appear complete 
to w 13 mag. The integration times were set so that 
the median polarization uncertainty, ap, would be less 
than 1% for stars brighter than H = 14. As shown in 
Figure [2l this goal was not quite achieved, with a me- 
dian polarization uncertainty of 1.3% at H — 14. The 
shape of the curve in Fig. [2] indicates that polarimetric 
uncertainties are limited by photon noise. 

The polarization efficiency {P/Ay) towards each star 
carries information about the polarization environ- 
ment and polarization mechanism (iGoodman et al.lll992l : 
iLazarian et"allll997l: iWhittet et al.ll2008D . Because these 
observations span from the Galactic midplane to 75 de- 
grees above the Galactic plane, a change in P/Ay could 
indicate a change in the polarization or grain alignment 
mechanism. The extinction towards each star can be es- 
timated from the o bserved (H-K) 2 MASS colors and the 
intrinsic colors from lBessel fc Bretti ([1988). Assuming all 
stars have an intrinsic (H-K) color of 0.1 ±0.1 mag, E(H- 
K) values can be converted to extinctions at V-band, Ay, 
as: 

Av - t'Eh-k, (1) 

assum ing r = 16 for the diffuse ISM (jWhittet et all 
119961 ). Considering only stars with measured 2MASS 
(H-K) colors, 2MASS photometric SNR > 3 for H and 
K, and polarimetric SNR > 3, the resulting photomet- 
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Table 2 

Starlight Polarization Observations 
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0.000 


99.990 


0.000 


99.990 



Note. — Table[2]is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. Stars with 
missing 2MASS photometry are listed as having apparent magnitudes of 0.00 mag and 
uncertainties of 99.99 mag 



ric and polarimetric extinction distributions are shown in 
Figure [31 The median extinction of the photometric and 
polarimetric stars are Ay = 2.6 and 1.8, respectively. In 
Fig. |3l the photometric and polarimetric distributions 
appear to have the same overall shape, indicating that 
they are drawn from the same parent distribution and 
that the polarimetric selection criteria does not intro- 
duce bias to that sample. The observed distribution of 
polarization with extinction was examined for trends in 
the polarization efficiency {P/Av) with Galactic latitude, 
but none were seen. This implies no detectable change 
in the polarization mechanism with Galactic latitude. 

In Figure HI the observed GPA of each significantly 
and marginally detected stellar polarization is shown, as 
a function of Galactic latitude. The typical GPA un- 
certainty is 5° for the significant detections and 21° for 
the marginal detections. Alignment with the Galactic 
plane would occur at GPA = 90°, as shown by the dot- 
ted line. Fairly good alignment is seen near the Galactic 
midplane (b ^ 0°), in agreement with the findings from 
optical polarimetry (Heiles 2000). To quantify the GPA 
for each field and how it varies with Galactic latitude, the 
GPA weighted means and uncertainties were calculated 
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Figure 3. Distributions of the estimated extinctions towards pho- 
tometric (dotted line) and significantly polarized (solid line) stars, 
as described in the text. 
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Figure 4. Measured Galactic position angle (GPA) as a func- 
tion of Galactic latitude for stars showing significant and marginal 
starlight polarizations. The dotted line at GPA = 90° represents 
perfect alignment of the magnetic field parallel to the Galactic 
plane. 

for each field and are shown in Figure [51 The error bars 
in the figure represent the weighted ±la uncertainties in 
the mean GPAs. 

tion CDF for stars exhibiting significant polarization 
[P/crp > 3) for all latitude bins is shown in Figure[Sl The 
1st, 2nd, and 3rd quartiles occur at 1.04%, 1.54%, and 
2.38%, respectively, and are shown by the dashed lines 
in Fig. m The smallest significantly detected polariza- 
tion was 0.19%, which is comparable to the polarization 
sensitivity limit of these observations. The polarization 
lower limit arises from both p hoton noise and calibration 
uncertainties, as discussed in [Clemens et alj (|2012bf ). 

3. COMPARISON WITH SIMULATED OBSERVATIONS 

Comparisons of these observations with predictions of 
the large-scale structure of the Galac tic magnetic field 
can constrain the magnetic geometry. iPavell (|2011t ) sim- 
ulated all-sky NIR starlight polarization observations for 
several magnetic field geometries. These included: three 
AO (disk antisymmetric), an d three SO (disk symmet - 
ric) magnetic field geom etries (iFerriere fc Schmitt]l2000[ l: 
three DEHO geometries (jMoss et al.ll2010t ): and three an- 
alytic axisymmetric magnetic fields with different pitch 
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Figure 5. Weighted mean GPA for each observed field, versus 
Galactic latitude. The error bars represent the weighted uncer- 
tainties of the means for each field, and the dotted line represents 
perfect alignment parallel to the Galactic plane. The dashed box 
identifies the Galactic latitude range that is most diagnostic for 
constraining the Galactic magnetic pitch angle (see text in §4.2). 
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Figure 6. The observed starlight polarization percentage CDF for 
all 386 high-significance {P/ap > 3) stars, drawn from all latitude 
bins. The dashed lines show the locations of the 1st, 2nd, and 3rd 
quartiles of the significantly detected polarization values. 

angles (a = 0°, —11.5°, and —24°). In a trailing spi- 
ral galaxy, since the magnetic field is tied to the gas 
and not to the spiral density pattern, differental rota- 
tion c auses the and B^, compon ents to have opposite 
signs (jKrasheninnikova et al.l [19891) . The pitch angle is 
defined hy p = ta n~^(Bj. / Bj,), a nd therefore the pitch 
angle is negative ()Beck et al.l [l996) . SO and AO mag- 
netic fields are two common axisymmetric predictions 
for the large-scale structure of the Galacti c magnetic 
field. DEHO magnetic fields, pr edicted by iMoss et al.l 
(|2010( ) and possibly observed by iSun et al.l (|2008D . in- 
voke a Galactic wind that connects independent halo and 
disk dynamos. Analytic axisymmetric magnetic fields 
(consisting o f only B^ and Bj. components) were used by 
iPavell (|2011| ) to predict the effects of a spiral-based Galac- 
tic magnetic pitch angle on the observed polarization dis- 
tribution, and do not represent a magnetic field geome- 
try with a known physical basis. All of these magnetic 
field models were com bined with two different empiri- 
cal dust distributions (ISpergel. Malhotra. fc Blitd [19961 : 



iDrimmel fc Spergell I2001D vielding a total of 24 sets of 
predictions, as summarized in the first three columns of 
Table m In addition to an all-sky map for each predic- 
tion, an £ = 150° GPA with h plot was shown, all-sky 



polarization CDFs were plotted, and the 1st, 2nd, and 
3rd all-sky CDF quartile values were tabulated for syn- 
t hetic stellar s amples complete io H — lA. 

IPavell (j2011l) suggested two tests were most diagnos- 
tic for comparing the predictions with observations: (1) 
comparing the CDFs of the observed and predicted 
starlight polarizations, and (2) comparing the change 
in GPA with Galactic latitude between the observations 
and predictions. Each test is applied and discussed be- 
low. 

3.1. Cumulative Distribution Functions 



The predicted CDFs in IPavell (|201lD were normalized 
in that work by the maximum predicted polarization, be- 
cause of uncertainties in the degree of polarization caused 
by the unknown magnetic alignment and dust polariza- 
tion efficiencies. If these factors do not change signifi- 
cantly along or between different lines of sight, then the 
predicted GPA is unaffected and the normalized CDFs 
are related to the actual polarization CDFs by a multi- 
plicati ve fact o r. 

The IPavell (|2011| ) predictions were generated for the 
whole sky, while the observations presented here cover 
only a small region of sky. Since the predicted degree of 
polarization strongly depends on Galactic latitude and 
longitude, a subsample was drawn from the all-sky pre- 
dictions which included only those directions correspond- 
ing to observations reported here. To account for the ob- 
servations at b = ±1.25°, ±2.5°, and —6° that have no 
corresponding |Pavel| (|2011l ) predictions, new simulations 
towards these direct ions were made, following the proce- 
dure in IPavell (|2011[ ). and included in the subsample. 

While the shapes of the variation of GP A with Galac- 
tic latitude for all 24 models presented in IPavell (j2011[ l 
are unchanged, the new subsample predicted polariza- 
tion CDFs are, a s wou ld be expected, different from 
the all-sky 'PaveH (|201lD CDF predictions. In Figure H 
the new subsample CDF predictions (solid lines), con- 
taining 332 stars, are shown along wit h mod i fied a ll-sky 
CDF predictions (dotted fines) from IPavell (|2011[ ). In 
the original all-sky predictions, a population of distant, 
high polarization stars exists towards the Galactic cen- 
ter (1^1 < 20°, 6 — 0°) that are not characteristic of the 
rest of the sk y. Th e presence of these stars was identi- 
fied in lPa^ |2nTl as a discontinuity in the CDFs pre- 
sented there. These very high polarization stars were re- 
moved from the all-sky CDFs shown in Fig. [7| (as dotted 
lines) . Significant differences remain between the shapes 
of the modified original and subsample CDFs. The in- 
clusion of four additional sightlines near the Galactic 
plane (6 = ±1.25°, ±2.5°) to the subsamples preferen- 
tially adds higher polarization stars, because of the large 
dust columns expected in the plane, which increases the 
cumulative probability at higher polarizations. This re- 
sults in a lower cumulative probability at small P/Pmax 
and higher probability at large P/Pmax, seen in all of the 
CDFs as a general flattening of the subsample curves. 

The multiplicative scale factors between each of the 
new predicted subsample CDFs and the observed CDF 
must also be calculated. Simple comparison of normal- 
ized predicted subsample CDFs to the observed CDF is 
unreliable because the maximum observed polarization 
value (used to normalize all other values) was based on a 
single star in the high end tail of the polarization prob- 
ability distribution. A single extreme polarization value 
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Table 3 

Model Parameters 



Model 


M^Sygnetic 


Dust 


K-S Test 




Fittin 


g ConstcintrS*^ 




AGFA 


Nuinber 


Field 


Model^ 


Probability 


A 


B 


C[xlO-''] 


D[xl0"5] 


Ides] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


SO reference run'' 


DS2001 


0.195 


90.1 


-0.72 


-1.53 


4.00 


7.14 ±0.94 


2 


SO reference run'' 


SMB96 


0.079 


90.1 


-0.71 


-0.98 


3.95 


6.99 ± 0.94 


3 


SO reference with alpha quenching' 


0.252 


DS2001 


90.0 


-0.58 


2.19 


3.35 


3.69 ± 0.94 


4 


SO reference with alpha quenching' 


0.094 


SMB96 


90.0 


-0.58 


2.67 


3.36 


3.60 ± 0.94 


5 


SO with vacuum BC' 


DS2001 


0.270 


90.0 


-0.57 


1.60 


3.22 


4.54 ±0.94 


6 


SO with vacuum BC' 


SMB96 


0.105 


90.0 


-0.57 


2.08 


3.24 


4.46 ± 0.94 


7 


AO reference run^ 


DS2001 


0.241 


-0.20 


0.043 


43.4 


9.76 


78.12 ± 1.03 


8 


l^\J 1 CiCl Cil^C 1 LIU 


SMB9fi 


0.055 


-0.03 


0.060 


13.0 


7.63 


77 r)4 -L 1 Kfi 


y 


AO reference with alpha quenching' 


JJozUUl 


U.Z41 


-U.z4 


U.Uoy 


41.0 


y.Do 


'7'7 VO J- 1 HQ 


10 


AO reference with alpha quenchmg" 


SMB9d 


0.084 


-0.12 


0.043 


4.19 


6.74 


75.91 ± 1.71 


11 


AO with vacuum BC' 


DS2001 


0.241 


-0.27 


0.047 


46.0 


10.14 


77.75 ± 1.08 


12 


AO with vacuum BC' 


SMB96 


0.063 


-0.19 


0.061 


12.6 


7.69 


75.94 ± 1.71 


13 


DEHO C„i„d = 0, Rc,,haio = 300" 


DS2001 


7.42 X lO"'' 


90.1 


-0.80 


-14.2 


3.31 


8.66 ± 1.15 


14 


DEHO C^ind = 0, Ra.halo = SOO'^ 


SMB96 


2.70 X 10"^ 


90.1 


-0.80 


-14.7 


3.24 


8.68 ± 1.17 


15 


DEHO Cyjind = 100, Ra,haio = 300" 


DS2001 


8.91 X 10"^ 


90.1 


-0.77 


3.87 


5.36 


8.69 ± 1.37 


16 


DEHO Cy^ind = 100, -R(j,ho!o = SOO'^ 


SMB96 


4.31 X 10-5 


90.1 


-0.78 


3.11 


5.35 


8.66 ± 1.39 


17 


DEHO C^i-ad = 200, Rc,haio = 300" 


DS2001 


7.25 X 10-5 


89.5 


-1.13 


-61.1 


0.64 


14.37 ± 1.28 


18 


DEHO Cy^ind = 200, Ret, halo = SOO'^ 


SMB96 


3.32 X 10-5 


89.5 


-1.13 


-61.8 


0.56 


14.34 ± 1.30 


19 


Ring, = 0" 


DS2001 


0.270 


90.0 


-0.54 


0.71 


2.95 


2.56 ± 0.93 


20 


Ring, = 0° 


SMB96 


0.167 


90.0 


-0.54 


1.23 


2.95 


2.49 ± 0.94 


21 


Ring, = 11.5° 


DS2001 


0.265 


90.0 


-0.31 


4.86 


1.89 


-3.58 ±0.93 


22 


Ring, = 11.5° 


SMB96 


0.222 


90.0 


-0.30 


5.25 


1.89 


-3.64 ±0.94 


23 


Ring, = 24° 


DS2001 


0.270 


90.0 


-0.080 


6.01 


0.86 


-9.66 ±0.93 


24 


Ring, = 24° 


SMB96 


0.156 


90.0 


-0.076 


6.34 


0.86 


-9.72 ±0.94 



DS2001 = 'Drimmel fc Spergell 120011 '): SMB96 = ISpergel. Malhotra, fc Blit2l HTMB) . 
' Ulrricrc & Schmitt ( 2Q00) 
" [Moss et al., (,2010) 

<i The fit is of the form CPA = A + Bxb + Cxb'^ + Dxb^. 

would affect all of the normalized starlight polarizations 
and change the cumulative distribution function. 

For each model comparison, the null hypothesis is that 
the subsample-predicted and NIR-observed polarization 
CDFs are drawn from the same parent distribution. If 
the probability that they are drawn from the same par- 
ent distribution is small, the null hypothesis can be re- 
jected. The Kolmogorov-Smirnov (K-S) test was used 
to determine both the appropriate scaling factor and to 
perform the test of the null hypothesis. By scaling the 
predicted subsample CDFs so that the maximum abso- 
lute difference between the observed and predicted CDFs 
was minimized, the K-S probability that the CDFs are 
drawn from the same parent distribution is maximized. 
This introduces a heavily conservative bias that overesti- 
mates the likelihood of agreement. If this overestimated 
probability is still significantly small, then that particu- 
lar model can be rejected. In this case, the overestimated 
probability is a strong, conservativly biased criteria for 
rejecting the null hypothesis. 

For sky locations matching the observations, 866 sim- 
ulated stars with H<14 mag make up the predicted sub- 
sample CDFs. The 332 observed stars with II<14 mag 
and significant polarizations make up the NIR-observed 
CDF. Many of the simulated stars, however, would ac- 
tually have polarization values below the NIR observa- 
tional limit and should be culled. To calculate the stars 
to ignore, the following procedure was used for each of 
the 24 subsample predictions. First, the entire predicted 
subsample polarization CDF was scaled so that the max- 
imum deviation between it and the observed polarization 
CDF was minimized. Next, model stars with scaled po- 



larizations below the observational polarization thresh- 
old (0.19%) were identified and removed. Then, a new 
predicted polarization CDF was calculated and scaled to 
the observations again. The resulting (overestimated) 
K-S probabilities for each model-observations compari- 
son are listed in column (4) of Table [31 These represent 
the probabilities that the best-scaled predicted CDFs 
were drawn from the same population as the observed 
CDF. The largest probability of all the models is only 
0.27 and is too low to identify any obvious good match. 
The smallest probabilities (strongest model rejections) 
are approximately 10^^ for the DEHO models, though 
further consideration is necessary before rejection and 
will be discussed below. 

3.2. GPA vs. Galactic Latitude 

The second observational test proposed by iPavell 
was based on the dependence of the polarization 
CPA on Galactic latitude. The observational data for 
this test were shown in Fig. Weighted means and un- 
certainties in the mean for each latitude bin were shown 
in Fig. [51 To facilitate comparison with the observations, 
the runs of predicted CPAs, for each model at £ = 150°, 
as a function of Galactic latitude (see Figs. 4-7 in IPavell 
120111 ) were fit by a third order polynomial: 

GPA = A + Bxb + Cxb^ + Dxb^ (2) 

where GPA is the predicted Galactic polarization posi- 
tion angle; b is the Galactic latitude of the prediction; 
and A, B, C, and D are constants to be fit. The parame- 
ters from these fits are listed in columns (5)- (8) in Table 
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Figure 7. Predicted starliglit polarization CDFs for tlie subse t of di r ection s corresponding to the observations reported here (solid lines). 
These curves are different from the all-sky curves presented in lPavell I I2011I , dotted lines; modified as described in the text), especially at 
low and high polarizations because of the limited sampling of the sky by the new observations and the latitude and longitude dependence 
of the polarizations. 

[31 The choice of a third-order polynomial was based on 
an F-test. 

The observed GPAs and two example model polyno- 
mial fits are plotted in Figure [H Model f (top panel) 
is an example of a fair match between the predictions 
and observations and Model 7 (bottom panel) is an ex- 
ample of a poor match. For each latitude bin having 
observations, the difference between the predicted mean 
GPA and observed mean GPA was calculated. These are 
shown in Figure |9] for all models. The weighted obser- 
vational mean GPA uncertainties and the predicted sub- 
sample GPA dispersions are combined in quadrature to 
estimate the uncertainties for these differences. In Fig. 
[21 Models 1-6 and 13-24 show fair agreement between 
predictions and observations. Models 7-12 (AO magnetic 
field geometries), however, do not agree with the obser- 
vations, with most latitude bins having discrepancies of 
AGFA > 80°, essentially orthogonal to the predicted 
GPAs, and therefore are not shown. 

To reject specific models, it is necessary to quantify 
these differences. For each model, a weighted average 
of the differences between the model's predicted GPA 
and observed GPA in each latitude bin, as shown in Fig. 
[9l was calculated and the weighted uncertainty in this 
average difference was obtained by propagation. These 
values and uncertainties are listed in column (9) in Table 
[3 While ah SO, AO, and DEHO models can be rejected 
at the 3(7 level, the AO models are particularly discrepant 
(>40ct) since they predict polarizations essentially per- 
pendicular to the Galactic plane. 

4. DISCUSSION 




-80 -60 -40 -20 20 
Galactic Latitude [deg] 

Figure 8. Two example model polynomial fits (dashed lines) plot- 
ted over the observed GPAs (points) as a function of Galactic lati- 
tude. For clarity, only the significantly detected (P/crp > 3) polar- 
izations are plotted, and are identical in both panels. The SO-type 
Model 1 (top) is an example of fair agreement between the pre- 
dictions and observations. The AO- type Model 7 (bottom) is an 
example of poor agreement, with an approximately 90° difference 
between the predictions and observations. 

4.1. Large-Scale Magnetic Field Symmetry 

Using the two tests proposed bv iPavell (12011 1 and the 
observations reported here, all of the magnetic field mod- 
els taken from the literature can be formally rejected as 
deviating by more than 3a from the observations. The K- 
S analysis of CDF values was inconclusive about rejection 
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Figure 9. The differences between the observed mean GPAs and predicted subsample mean GPAs in each latitude bin for Models 1-6 
and 13-24. Models 7-12 are not shown because the differences for all latitude bins are greater than 60°. The error bars show the itlo" 
uncertainties for each difference. The effect of magnetic pitch angle on AG PA for |b| > 10° is seen in Models 19-24 where the pitch angle 
(listed to the right of each plot) is varied. 



of the predicted polarization values. None of the proba- 
bilities was high enough to suggest one set of models was 
the best choice, and (except for the DEHO models) the 
probabilities were not small enough to definitively reject 
classes of models outright. The differences between the 
predicted and observed polarization GPAs as a function 
of Galactic latitude are significant enough to reject all of 
the magnetic field models. However, the full interpreta- 
tion of these results requires a more careful analysis. 

The initial results should not be completely surpris- 
ing, since the actual Galactic magnetic field is unlikely 
to be as simply represented as the simulated magnetic 
field models. The dynamo-driven magnetic field m od- 
els have a resolution of ^^200 pc ([Moss et all l2010f) or 
^400 pc (jFerriere fc SchmittI 120001 ) . so anv turbulence 
or magnetic disturbances on smaller scales could not 
be included. The analytic axisymmetric models, with 
effectively infinite resolution, also lack turbulence. If 
turbulence creates random magnetic field components, 
then along a line of sight the starlight is expected to 
suffer some depolarizatio n as it passe s through multiple 
rando mly oriented cells ponesi [l989l ) . lOhno fc Shibatal 
(|1993l ) calculated a coherent magnetic cell size of 10 to 



100 pc using Galactic pulsar rotation and dispersion mea- 
sures, implying that starlight fro m dista n t star s will gen- 
erally pass through many cells. iPavell ()2011l) also dis- 
cussed uncertainties in the predicted degree of polariza- 
tion arising from changes in the polarization efficiency, 
though they are less likely to affect the GPAs. This, 
combined with uncertainties from the calculated scale 
factors, creates uncertainty in the absolute probability 
returned by the K-S test, though the relative probabili- 
ties may still provide insight into ranked likelihoods. 

Considering the relative K-S probabilities in Table [3l 
the SO, AO, and analy t ic axi symmetric models using 
the iDrimmel fc Spergel (|2001[ ) dust distributions have 
the highest probability (> 0.20) of being drawn from 
the same parent distribution as the observations. The 
DEHO models all produce polarization CDFs that are 
poor matches to the observed polarization CDF. The 
choice of dust model causes a significant difference in 
the probabilities for all models. T he SO, AO, and ana- 
lytic a xisymmetric models using the IDrimmel fc Spergell 
((2OOII ) dust distribution have larger (i.e., more similar 
to t he observations) probabilities tha n the models using 
the iSpergel. Malhotra. fc Blitd Hl996f ) dust distribution. 



Testing Galactic Magnetic Fields 



9 



This i s not surprising given that the IDrimmel fc Spergell 
(|2001h dust model accounts for more details (e.g., spiral 
arr ns, local arm, Galactic wa r p) tha n the axisymmet- 
ric iSpergel. Malhotra. fc Blitg (|1996D dust model. The 
DEHO models, as shown by their K-S probabilities, are 
unable to reproduce the observed polarizations. 

The amount of polarization increases as starlight 
passes through the regular (uniform) magnetic field, but 
depolarizes as it passes through cells of randomly ori- 
ented magnetic field. The interplay between regular and 
random magnetic fields is strongly dependent on the na- 
ture of turbulence, and beyond the scope of this pa- 
per. Because of the lack of depolarization in the simu- 
lations, all simulated stars will show larger polarizations 
than would actual stars. The exact effect on the sim- 
ulated degrees of polarization requires a detailed study 
of the ratio of power in the random and regular mag- 
netic fields in the diffuse ISM, analogous to the studies 
that ha ve been conducted in higher density molecular 
clouds (iFalceta-Goncalves et al.l 120081 : IHildebrand et al.l 
I2009t IHoude et all 120091 ). The changes to the degree 
of polarization caused by depolarization may be large 
enough to change the predicted polarization CDF and 
account for the significant differences measured by the 
K-S test. 

Given the uncertainties associated with the predicted 
degree of starlight polarization, the distribution of GPA 
vs. Galactic latitude may serve as a better test for con- 
straining the large-scale structure of the Galactic mag- 
netic field. 

The AO (disk-odd) Galactic magnetic fields, models 7- 
12, show huge differences between the observed and pre- 
dicted GPAs towards i = 150° and are therefore not 
shown in Fig. [S] In typical AO magnetic field models, 
the poloidal magnetic field has a dipolar configuration 
and the toroidal magnetic field consists of antisymmet- 
ric tori of magnetic fiux above and below the Galac- 
tic disk with opposite field directions. To sustain this 
structure, there must be a toroidal null point approx- 
imately in the Galactic plane where the toroidal mag- 
netic field goes to zero. At this point, only the poloidal 
magnetic field is expressed, which is orientated roughly 
in the Galactic North-South direction, orthogonal to the 
Galactic plane as seen in the simulations. Also, most of 
the toroidal magnetic flux in the specific AO magnetic 
fie ld geometries used for the se simulations (e.g.. Fig. 10 
in iFerriere &: Schmit^ |2000i for Model 7) is contained 
within the solar circle and would not be probed when 
looking toward the outer Galaxy. Based solely on the 
GPA discrepancies, these AO models can be rejected. 

The rejection of AO magnetic field geometries is at 
odds wit h previous work utilizing rotation measure stud- 
ies (e .g., Andreasvan fc Makarovl 119881: iHan et al.l 119971 
l2003t [Sun et al. 2008), which probe the Galactic mag- 
netic field over kiloparsec length lines of sight. Evidence 
for an antisymmetric magnetic field is seen in antisym- 
metric rotation measures above and below the Galac- 
tic plane towards the Galactic center. This implies that 
the dominant toroidal magnetic field dire ction is reversed 
above and below the Galactic plane. iWolleben et al.l 
(|2010| ). however, has recently shown that this observed 
asymmetry may be dominated by a nearby (100 pc), 
northern-sky, HI bubble and that the observed asymme- 
try may not reflect the actual large-scale Galactic mag- 
netic field. By using a different method (starlight po- 
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Figure 10. The average difference between the observed and pre- 
dicted GPA distributions at ^ = 150°, -60° < 6 < -15° as a 
function of magnetic field spiral pitch angle for the analytic ax- 
isymmetric model. The horizontal dashed line shows where perfect 
agreement between the observations and predictions occurs. The 
error bars represent the itlo" uncertainties on the mean differences. 
Only models with error bars enclosing APA = 0° are shown. The 
diagonal solid lines trace the ztlcr errors of the tested models. The 
vertical solid line is the best fit magnetic pitch angle and the grey 
region shows the estimated uncertainty on that value, as described 
in the text. 

larimetry instead of Faraday rotation) away from known 
magnetic disturbances, the rejection of the AO magnetic 
field is robust. 

The analytic axisymmetric models (Models 19-24) 
show the best overall matches, both in degree of polar- 
ization and GPA. The difference between the observed 
and predicted GPAs suggests that —24° is too large to 
match the observations. Based on the quality of these 
matches, a magnetic pitch angle between 0° and —11.5° 
provides the best fit to the data. 

4.2. Magnetic Pitch Angle 

To place tighter constraints on the magnetic pitch an- 
gle, additional analytic axisymmetric model polarization 
subsample predictions were generated for pitch angles of 
9 — 0° to —24° in steps of 0.5°. The average discrepan- 
cies between the observed and predicted GPAs were cal- 
culated for the data in the zone —60° < b < —15°, shown 
by the dashed box in Fig. [5] (all latitudes were used in 
the earlier analysis). Latitude bins close to the Galac- 
tic plane were not used for two non-independent reasons: 
the average predicted GPAs for |6| < 10° are not signif- 
icantly affected by the changing pitch angle (as seen in 
the rightmost column of Fig. ^ , and the uncertainty on 
the GPA differences near the midplane are much smaller 
than those at larger latitudes. For a properly weighted 
mean, the change in uncertainty with latitude will cause 
the non-diagnostic, low latitude bins to dominate the av- 
erage while the diagnostic, large latitude bins carry little 
weight. Also, the observations at 5 = —65° and —75° 
were not used because their GPAs were not consistent 
with any of the models and include only one significant 
starlight polarization measurement (see Table [T]). 

The average difference between the predicted and ob- 
served GPAs, and ±lcr uncertainties, as a function of 
magnetic pitch angle, are shown in Figure lTOl In this Fig- 
ure, the six models whose error bars enclose AGFA — 0° 
(horizontal dotted line) are shown. The models not 
shown follow the same trend away from AGFA = 0° with 
larger and smaller model magnetic pitch angles. The di- 
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agonal solid lines trace the ilcr error bars for all of the 
models. The best fit magnetic pitch angle is between 
9 = —5.5° and —6.0°, shown by the vertical solid line in 
Fig. [TUl To estimate the uncertainty on this magnetic 
pitch angle, the upper and lower la uncertainty limits 
for all analytic axisymmetric models were each fit with 
second order polynomials (shown as dotted lines in Fig. 
[T0|) . The pitch angles enclosed by these fits are shown by 
the grey region in Fig. [101 Based on these new observa- 
tions and predictions, the Galactic magnetic pitch angle 
is 6* = -6 ± 2° towards i = 150°. 

Previous estimates of the Galactic mag netic pitch an - 
gle have retu rned similar result s: 
-8.2 ° ± 0.5° (IHan fc Qiaol lTool . 

[1991), and -8° ± 1° (|Be'cni200l . Unlike these previous 



estimates, which typically use data from large areas of 
the sky, this result was obtained from one Galactic lon- 
gitude. Additional measurements in the outer Galaxy 
may show pitch angle variations with Galactic longitude 
that are missed when averaging over large areas. 

5. CONCLUSIONS 

New NIR starlight polarimetry observations were ob- 
tained toward the outer Galaxy that trace the orienta- 
tion of the Galactic magnetic fiel d. Th e se obs ervations, 
combined with predictions from iPavell ()201lD and new 
simulations presented here, constrain important aspects 
of the geometry of the large-scale Galactic magnetic field. 
The Kolmogorov-SmirnofF (K-S) test probabilities of the 
cumulative distribution functions (CDFs) of the polariza- 
tion percentages and the Galactic position angle (GPA) 
differences can be used to formally (Str) reject all of the 
models tested here, though relative rankings may also be 
useful. Key results from this study include: 

1. The AO model predictions are strongly inconsis- 
tent with the observed distribution of NIR starlight 
Galactic position angles (GPAs) and are thereby 
rejected. W hile here limited to a set o f sample AO 
models from iFerriere fc Schmitti ()2000D . the ability 
of any AO-type magnetic field model to reproduce 
the observed GPA distribution in the outer Galaxy 
is doubtful. 

2. The GPA distributions of the SO, DEHO, and ana- 
lytic axisymmetric models typically differ from the 
NIR starlight observations by only a few degrees. 
Since the model magnetic fields do not include any 
turbulence or systematic motions, such differences 
are to be expected. The general agreement shows 
that disk-even (SO, DEHO, and analytic axisym- 
metric) magnetic fields better represent the Galac- 
tic magnetic field than AO models. 

3. Based on the CDFs of observed and predicted de- 
grees of starlight polarization and their K-S prob- 
abilities (see Table [3|), the Dis k-Even, Halo-Odd 
(DEHO) magnetic field models (jMoss et al.ll2010( ) 
are especially discrepant. This result, and all of the 
K-S probabilities, may be affected by uncertainties 
in predicting the degree of starlight polarization. 

4. Even if the predicted DEHO polarization CDF was 
not discrepant, the observed distribution of GPAs 
would still be inconsistent, though at a similar level 
as for the SO and analytic axisymmetric models. 



5. New simulations presented here, with different 
spiral-type magnetic pitch angles, when matched to 
the new observations, constrain the Galactic mag- 
netic pitch angle to —6° ± 2°. This estimate is 
for a single Galactic longitude and agrees well with 
previous estimates of optical starlight polarization 
jHcilcs 1996) and radio Fara day ro tation (jValled 
[1988.; Han fc Qiaol[T99l lBecS[2007l ) that each av- 
erage over large areas of the sky. 

This work highlights the utility of NIR polarimetry 
and its ability to place strong constraints on the Galactic 
magnetic field. Future application of NIR polarimetry to 
larger areas of the sky will provide additional constraints 
for Galactic dynamo theory on large scales and star for- 
mation on small scales. 
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